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ABSTRACT 

We present here our first series in studying the astrophysical parameters of the 
open cluster "Ruprecht 15" using PPMXLL 1 ! database. In this context, the pho- 
tometric, astrometry and statistical parameters for this cluster (limited radius, 
core and tidal radii, distances, membership, reddening, age, luminosity function, 
mass function, total mass, and the dynamical relaxation time) are determined 
for the first time. 

Subject headings: open clusters and associations - individual: Ruprecht 15 - 
astrometry - Stars - astronomical databases: catalogues. 

1. Introduction 

Open star clusters (OCs) play an important role in studying the formation and evolution 
of the Galactic disk and the stellar evolution as well. Also, OCs are the most suitable 
objects for studying the space-age structure of the Galactic disk. The fundamental physical 
parameters of OCs, e.g. distance, reddening, age, and metallicity are necessary for studying 
the Galactic disk. The Galactic, radial and vertical, abundance gradient can be studied by 
OCs (Hou et al. 2000; Chen et al. 2003; Kim and Sung 2003; Tadross 2003; Kim et al. 2005). 
They are excellent probes of the Galactic disk structure (Janes & Phelps, 1994; Bica et al. 
2006). The strong interest of OCs results come from their fundamental properties. Among 
the 1787 currently OCs, most than half of them have been poorly studied or even unstudied 
up to now, Piatti et al. (2011). The current paper is thus part of our continuation series 
whose goal is to obtain the main astrophysical properties of previously unstudied OCs using 
modern databases. The most important thing for using PPMXL database lies in containing 
the positions, proper motions of USNO-Bl.Oj and the Near Infrared (NIR) photometry of 
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the Two Micron All Sky Survey (2MASSj3, which let it be the powerful detection of the star 
clusters behind the hydrogen thick clouds those concentrate on the Galactic plane. 

The only available information about Ruprecht 15 (hereafter Ru 15) are the coordinates 
and the optical apparent diameter, which were obtained from WEBDA0 site and the last 
updated version of DIAS0 collection (version 3.0, 2010 April 30). This cluster is situated in 
the Southern Milky Way at 2000.0 coordinates a = 07 h 19 m 34 s , 5 = -19° 37' 30", £ = 
233.54°, b = —2.896°, and its diameter is about 2.0 arcmin. Fig. 1 represents the blue image 
of Ru 15 as taken from Digitized Sky Survey (left panel), while the right panel represents 
J-image of the cluster as taken from Interactive 2MASS Image Servicqj. 

This paper is organized as follows. PPMXL data extraction and preparation are pre- 
sented in Section 2, while the data analysis and parameters estimations are described in 
Sections 3. Finally, the conclusion is devoted to Section 4. 



2. PPMXL Data Extraction and preparation 

The astrophysical parameters of the investigated cluster have been estimated using the 
PPMXL Catalogue of Roeser et al. (2010). It is combining the proper motion of USNO-B1.0 
and NIR JHK S pass-band of 2MASS databases. 

USNO-B1.0 of Monet et al. (2003) is a spatially unlimited catalogue that presents 
positions, proper motions, and magnitudes in various optical pass-bands. The data were 
obtained from scans of Schmidt plates taken for the various sky surveys during the last 
50 years. USNO-B is believed to provide all-sky coverage, completeness down to V = 21 
mag. It is noted that, based on the proper motion measurements, stars with large proper 
motions are likely to be foreground stars instead of cluster members. Background stars 
cannot readily be distinguished from members by proper motions. Nonetheless, identifying 
foreground stars is useful in cleaning up the colour-magnitude diagrams and determining 
the field star contamination. So, USNO-B is a very useful catalogue, which gives us an 
opportunity to distinguish between the members and background/field stars. 

On the other hand, the NIR 2MASS photometry provides some direct answers to ques- 
tions on the large scale structure of the Milky Way Galaxy and the local Universe. It provides 
J, H, and K s band photometry for millions of galaxies and nearly a half-billion stars (Carpen- 
ter, 2001). 2MASS observations are obtained using two highly automated 1.3-m telescopes, 
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one at Mount Hopkins in Arizona (Northern Survey) and the other at Cerro Tololo in Chile 
(Southern Survey). Each telescope was equipped with three-channel camera, each channel 
consisting of 256 x 256 array HgCdTe detectors. It is uniformly scanning the entire sky in 
the three NIR bands J (1.25/im), H (1.65/iin) and K s (2.17//m). This survey has proven to 
be a powerful tool in the analysis of the structure and stellar content of open clusters (Bica 
et al. 2003, Bonatto & Bica 2003). The photometric uncertainty of the 2MASS data is less 
than 0.155 at K s ~ 16.5 magnitude which is the photometric completeness for stars with 
\b\ > 25°, Skrutskie et al. (2006). 

It is noted that Ru 15 is located near the Galactic plane (\b\ < 3°), therefore we expect 
significant foreground and background field contamination. The apparent diameter is less 
than 5 arcmin, hence the downloaded data has been extended to reach the field background 
stars, whereas the cluster dissolved there, so that we extracted the data to a radius of about 
10 arcmin. 

Comparing the data with Cutri et al.'s Point Source Catalogue (2003), we found that 
most stars' magnitudes have a "AAA" quality flag, which means the Signal Noise Ratio is 
SNR > 10, i.e. they have the highest quality measurements. Every star in our data has 3- 
colour photometric values J, H, K s mag; and proper motion (pm) values in right ascension (a) 
and declination (<5), i.e. (pm a cos 5 & pm 5) mas yr _1 . According to Roeser et al. (2010), the 
stars with proper motion uncertainties > 4.0 mas yr _1 have been removed. In this context, 
to get a net worksheet data for Ru 15, the photometric completeness limit has been applied 
on the photometric pass-band 2MASS data to avoid the over-sampling at the lower parts of 
the cluster's CMDs (cf. Bonatto et al. 2004). The stars with observational uncertainties > 
0.20 mag have been removed. Pm vector point diagram (VPD) with distribution histogram 
of 2 mas yr _1 bins for (pm a cos 5) and (pm S) have been constructed as shown in Fig. 2. 
The Gaussian function fit to the central bins provides the mean pm in both directions. All 
data lie at that mean ±1 o (where o is the standard deviation of the mean) can be considered 
as astrometric probable members. In addition, the stellar photometric membership criteria 
are adopted based on the location of the stars within ±0.1 mag around the ZAMS curve in 
the colour magnitude diagrams (CMDs). 



3. Data Analysis 

3.1. Cluster's Centre and Radial Density Profile 

The centre of any cluster can be roughly estimated by eye, but to determine the centre's 
coordinates of Ru 15 more precisely, we applied the star-count method to the whole area of 
the cluster. All the 10 arcmin data that obtained around the adopted centre are dividing 
into equal sized bins in a and 5. The cluster centre is define as the location of maximum 
stellar density of the cluster's area. The Gaussian curve-fitting is applied to the profiles of 
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Fig. 1. — Left panel represents the blue image of Ru 15 as taken from Digitized Sky Surveys, 
while the right panel represents J-image of the cluster as taken from Interactive 2MASS 
Image Service, the small open circle indicates the cluster's central part. North is up, east on 
the left. 




Fig. 2. — Proper motion vector point diagram VPD of Ru 15 after avoiding all data with 
pm errors > 4 mas yr _1 . Histograms of 2 mas yr" 1 bins in both directions are drawn. The 
Gaussian function fit to the central bins provides the mean pm a cos 5 = -0.84 ±0.1 mas 
yr _1 and pm 5 = 6.7 ± 0.09 mas yr _1 . 
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star counts in a & 5 respectively as shown in Fig. 3. The differences between our estimated 
centre and the obtained one of Webda are shown in the figure. 

To establish the radial density profile (RDP) of Ru 15, we counted the stars within 
concentric shells in equal incremental steps (r < 1) arcmin from the cluster centre. We 
repeated this process for 1 < r < 2 up to r < 10 arcmin, i.e. the stellar density is derived 
out to the preliminary radius of the cluster. The stars of the next steps should be subtracted 
from the later ones, so that we obtained only the amount of the stars within the relevant 
shell's area, not a cumulative count. Finally, we divided the star counts in each shell to the 
area of that shell those stars belong to. The density uncertainties in each shell was calculated 
using Poisson noise statistics. Fig. 4 shows the RDP from the new centre of Ru 15 to the 
maximum angular separation of 5 arcmin where the stability of density has been reached. To 
determine the structural parameters of the cluster more precisely, we applied the empirical 
King model (1966). The King model parameterizes the density function p(r) as: 



where fb g , fo and r c are background, central star density and the core radius of the cluster 
respectively. From the concentration parameter c, defined as c = (Rum/R C ore), Nilakshi et 
al. (2002) concluded that the angular size of the coronal region is about 6 times the core 
radius. Maciejewski & Niedzielski (2007) reported that Ru m may vary for individual clusters 
between about 2R core and 7R core . In our case, we can see that Ru m =6.9R core , i.e. it lies 
within the previous vales. The cluster limited radius can be defined at that radius which 
covers the entire cluster area and reaches enough stability with the background density, 
i.e. the difference between the observed density profile and the background one is almost 
equal zero. It is noted that the determination of a cluster radius is made by the spatial 
coverage and uniformity of PPMXL photometry which allows one to obtain reliable data on 
the projected distribution of stars for large extensions to the clusters' halos. On the other 
hand, the concentration parameter seems to be related to cluster age, i.e. for clusters younger 
than about 1 Gyr, it tends to increase with cluster age. Nilakshi et al. (2002) notes that 
the halos' sizes are smaller for older systems. Finally, we can infer that open clusters appear 
to be somewhat larger in the near-infrared than in the optical data, Sharma et al. (2006). 
Knowing the cluster's total mass (Sec. 3.3), the tidal radius can be given by applying the 
equation of Jeffries et al. (2001): 



R t = 1.46 (M c ) 1/3 = 10.7 pc (2) 
where R t and M c are the tidal radius and total mass of the cluster respectively. 
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Fig. 3. — Estimating the cluster centre coordinates of Ru 15. The Gaussian fit provides the 
coordinates of highest density areas in a & 5 as 07 h 19 m 38 s & - 19° 37' 35" respectively. 
The differences between our estimated centre and the obtained one of Webda are 4 s in a 
and 5 in 5. 




Fig. 4. — The radial density distribution for stars in the field of Ru 15. The density shows a 
maximum at the centre p = 63 stars/arcmin 2 and then decreases down to p = 3 stars/arcmin 2 
at 2.2 arcmin, where the decrease becomes asymptotical at that point. The curved solid line 
represents the fitting of King (1966) model. Errors bars are determined from sampling 
statistics (1/y/N where N is the number of stars used in the density estimation at that 
point). The dashed line represents the background field density, where ft, g = 3.0 stars per 
arcmin 2 . The core radius r c = 0.32 arc min. 
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3.2. Colour-Magnitude Diagrams 

The main photometrical parameters (age, reddening and distance) are determined for 
our cluster by fitting the 2MASS ZAMS solar metallicity isochrones of Marigo et al. (2008); 
selected and downloaded from Padova isochronefl to the both CMDs (J, J-H & K s , J-K a ) 
of Ru 15. Several isochrones of different ages are applied; the best fit should be obtained 
at the same distance modulus for both diagrams, and the colour excesses should be obeyed 
Fiorucci & Munari (2003) 's relations for normal interstellar medium as shown in Fig. 5. 

The observed data has been corrected for interstellar reddening using the coefficients 
ratios = 0.276 and ^ = 0.176, which were derived from absorption rations in Schlegel 

et al. (1998), while the ratio ^ = 0.118 was derived from Dutra et al. (2002). 

Fiorucci & Munari (2003) calculated the colour excess values for 2MASS photometric 
system. We ended up with the following results: = 0.309±0.130, §^ = 0.485±0.150, 

where Rv = = 3.1. Also, we can de-reddened the distance modulus using these 

formulae: E g J = 0.887, E * a = 0.322, then the distance of the cluster from the Sun (R Q ) 
can be calculated. Therefore, (m-M)j =11.90 ± 0.10 mag (~ 1845 ± 85 pc), E(J-H) = 0.20 
± 0.05 mag, and (J-K s ) = 0.31 ± 0.07 mag. 

Once the cluster's distance R Q is estimated, then the distance from the galactic centre 
(R g ) and the projected distances on the galactic plane from the Sun (A & Y Q ) and the 
distance from galactic plane (Zq) can be determined, see Table 1. Fig. 6 illustrates a 
sketch-chart of the geometric distances for a cluster in the Galaxy, obtained from Tadross 
(2000). 



3.3. Luminosity, Mass Functions and the total mass 

Since the central concentration of the cluster is relatively low, the determination of the 
membership using the stellar RDP is difficult. Even though we determine the membership 
within a circle of radius r = 2 r c arcmin. By doing that, we obtained a more precise 
main-sequence in the CMDs. All of these stars are found very close to the main-sequence 
(MS) curve merged with some contaminated stars. These MS stars are very important in 
determining the luminosity and mass functions. 

The luminosity function (LF) gives the number of stars per luminosity interval, or in 
other words, the number of stars in each magnitude bin of the cluster. It is used to study the 
properties of large groups or classes of objects such as the stars in clusters or the galaxies 
in the Local Group. In order to estimate the LF of Ru 15 we count the observed stars in 
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Fig. 5. — The net NIR CMDs of Ru 15 for stars lying closely to the fitted isochrones and 
after removing all contaminated field stars. Age = 500 Myr, distance modulus = 11.9 mag 
(~ 1845 pc), E(J-H) = 0.20 mag, and (J-K s ) = 0.31 mag. 




Fig. 6. — A sketch-chart of the geometric distances for a cluster in the Galaxy. Let (d & 
D) be the angular and linear diameter in arc min and pc respectively, (-R©) the distance 
from the Sun, (R Q ) the sun's distance from the galactic centre (8.5 kpc), (Z Q ) the distance 
from the galactic plane, (R g ) the distance from the galactic centre, and (X Q & Y & ) the 
horizontal projected rectangular distances of the cluster from the Sun on the galactic plane; 
while (£ & b) are the galactic longitude and latitude of the cluster. 
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terms of absolute magnitude after applying the distance modulus as shown in Fig. 7. The 
magnitude bin intervals are selected to include a reasonable number of stars in each bin and 
for the best possible statistics of the luminosity and mass functions. From LF, we can infer 
that the massive bright stars seem to be centrally concentrated more than the low masses 
and fainter ones (Montgomery et al. 1993); the total luminosity is found to be ~ —4.1 mag. 

As known the LF and mass function (MF) are correlated to each other according to the 
Mass-luminosity relation. The accurate determination of both of them (LF & MF) suffers 
from some problems e.g. the field contamination of the cluster members; the observed 
incompleteness at low-luminosity (or low-mass) stars; and mass segregation, which may 
affect even poorly populated, relatively young clusters (Scalo 1998). On the other hand, the 
properties and evolution of a star are closely related to its mass, so the determination of the 
initial mass function (IMF) is needed, that is an important diagnostic tool for studying large 
quantities of star clusters. IMF is an empirical relation that describes the mass distribution 
(a histogram of stellar masses) of a population of stars in terms of their theoretical initial 
mass (the mass they were formed with). The IMF is defined in terms of a power law as 
following: 



m K M ' a (3) 



where is the number of stars of mass interval (M:M+dM), and a is a dimensionless 
exponent. The IMF for massive stars (> 1 M Q ) has been studied and well established by 
Salpeter (1955), where a = 2.35. This form of Salpeter shows that the number of stars in 
each mass range decreases rapidly with increasing mass. It is noted that our investigated 
cluster Ru 15 has MF slope ranging around Salpeter's value as shown in Fig. 8. 

To estimate the total mass of Ru 15, the mass of each star has been estimated from a 
polynomial equation developed from the data of the solar metallicity isochrones (absolute 
magnitudes vs. actual masses) at the age of the cluster (500 Myr). The summation of 
multiplying the number of stars in each bin by the mean mass of that bin yields the total 
mass of the cluster, which is about 390 M . 



3.4. Dynamical state 

The relaxation time (T re i ax ) of a cluster is defined as the time in which the cluster needs 
from the very beginning to build itself and reach the stability state against the contraction 
and destruction forces, e.g. gas pressure, turbulence, rotation, and the magnetic field (cf. 
Tadross 2005). T re \ ax is depending mainly on the number of members and the cluster diam- 
eter. To describe the dynamical state of the cluster, the relaxation time can be calculated 
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Fig. 7. — Spatial distribution of luminosity function of Ru 15 in terms of the absolute 
magnitude Mj. The colour and magnitude filters cutoffs have been applied to the cluster 
(dark area) and the field (white area) respectively. 
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Fig. 8. — The mass function of Ru 15. The slope of the initial mass function IMF is found 
to be T = —2.37 ± 0.20; with correlation coefficient of 0.82. 



- 11 - 



in the form: 

N 

Trelax 77 T cross , (4) 

8 m W 

where T cross = D/ay denotes the crossing time, N is the total number of stars in the 
investigated region of diameter D, and ay is the velocity dispersion (Binney & Tremaine 
1998) with a typical value of 3 km s _1 (Binney & Merrifield 1987). Using the above formula 
we can estimate the dynamical relaxation time for R 15, and then the dynamical-evolution 
parameter r can be calculate for the cluster by: 

age 

t = 7 , (5) 

^relax 

If the cluster's age is found greater than its relaxation time, i.e. r 3> 1.0, then the cluster 
was dynamically relaxed, and vice versa. In our case, Ru 15 is indeed dynamically relaxed, 
where r > 100. 



4. Conclusions 

Our procedure analysis has been applied for estimating the astrophysical parameters of 
yet unstudied open cluster "Ruprecht 15". Hence, we found that this cluster has real stellar 
density profile, its stars are lying at the same absolute distance modulus, reddening range, 
and having IMF slope around the Salpeter's (1955) value. On the other hand, the age of 
this cluster is found to be greater than the relaxation time, which infers that it is indeed 
dynamically relaxed. All parameters of this investigated cluster are listed in Table I. 

It is worthy to mention that, this publication made use of WEBDA, DIAS catalogues, 
and the data products from PPMXL database of Roeser et al. (2010). 
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Table 1: The present results of Ruprecht 15. 



Parameter 


The present result 


Center 


a = 07H9 m 38 s 




S = -19°37'35" 


c 


6.9 (see Sec. 3.1) 


pm a cos 5 


-0.84 ± 0.10 mas yr" 1 


pm 5 


6.7 ± 0.09 mas yr" 1 


Age 


500 ± 60 Myr. 


Metal abundance 


0.019 


E(J-H) 


0.20 ± 0.05 mag. 


E(J-K S ) 


0.31 ± 0.07 mag. 


E(B - V) 


0.65 ± 0.05 mag. 


Ry 


3.1 (see Sec. 3.2) 


Distance Modulus 


11.90 ± 0.10 mag. 


Distance 


1845 ± 85 pc. 


Rlim 


2.2' (1.20 pc.) 


Membership 


265 stars 


So 


63 ± 2 stars/arcmin 2 


fbg 


3 stars/arcmin 2 


Rc 


0.32' ±0.04 (0.20 pc) 


Rt 


10.7 pc. 


Rg 


9.7 kpc. 


x 


1095 kpc. 


Y 


-1482 kpc. 


Z 


-93 pc. 


T 


> 100 (see Sec. 3.4) 


Total Luminosity 


~ —4.1 mag. 


IMF slope 


T = -2.37 ± 0.20 


Total mass 


- 390 M Q 


Relaxation time 


2.7 Myr 



